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ABSTRACT 

A method is developed to compute the chord length distribution along a line 
which intersects a cellular Universe. The cellular Universe is here modeled by 
the Poissonian Voronoi Tessellation (PVT) and by a non-Poissonian Voronoi 
Tessellation (NPVT). The distribution of the spheres is obtained from com- 
mon approximations used in modeling the volumes of Voronoi Diagrams. 

We give analytical formulas for the distributions of the lengths of chords 
in both the PVT and NPVT. The astrophysical applications are made to the 
real Eso Slice Project and to an artificial slice of galaxies which simulates the 
2dF Galaxy Redshift Survey. 
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1. INTRODUCTION 
The study of the spatial periodicity in the spatial distribution of quasars 



started with Tifft (1973 1980 1995). A recent analysis of Bell & McDiarmid 



( 2006 ) , in which 46 400 quasars were processed, quotes a periodicity near Az 



-0.7. The study of periodicity in the spatial distribution of galaxies started 



with Broadhurst et al. (1990), where the data from four distinct surveys at 
the north and south Galactic poles were processed. He found an apparent 
regularity in the galaxy distribution with a characteristic scale of 128 Mpc. 
Recently, Hartnett ( 2009b|a ) quotes peaks in the distribution of galaxies with 
a periodicity near Az =0.01. More precisely, he found a regular real space 
radial distance spacings of 31.7 Mpc, 73.4 Mpc, and 127 Mpc. On adopting a 
theoretical point of view, the periodicity is not easy to explain. A framework 
which explains the periodicity is given by the cellular universe in which the 
galaxies are situated on the faces of irregular polyhedrons. A reasonable 
model for the cellular universe is the Poissonian Voronoi Tessellation (PVT); 
another is the non-Poissonian Tessellation (NPVT). Some properties of the 
PVT can be deduced by introducing the averaged radius of a polyhedron, R. 
The astronomical counterpart is the averaged radius of the cosmic voids as 
given by the Sloan Digital Sky Survey (SDSS) R7, which is R = ^-Mpc, 
see Pan et al. (2012). The number of intercepted voids, n v , along a line will 
be n„ = ~, where L is the considered length and I the average chord. On 
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Fig. 1. The number of galaxies as a function of the distance for a pencil beam 
catalog of Q, — 60 deg 2 . The curve was calibrated on the data of the 2dF Galaxy 
Redshift Survey (2dFGRS) which has Q = 1500 deg 2 . Adapted from Fi gure 4 in 



Zaninetti (20101 



assuming I = ny = |i. The astronomical counterpart of the line is the 
pencil beam catalog, a cone characterized by a narrow solid angle, see |Szalayl 



et al. (1993). The number of galaxies intercepted on the faces of the PVT 



will follow the photometric rules as presented in Zaninetti (20101, and Figure 
[T] reports the theoretical number of galaxies for a pencil beam characterized 
by the solid angle fl of 60 deg 2 . According to the cellular structure of the 
local universe, the number of galaxies as a function of distance will follow a 
discontinuous rather than a continuous increase or decrease. 
We can therefore raise the following questions. 

• Can we find an analytical expression for the chord length distribution 
for lines which intersect many PVT or NPVT polyhedrons? 

• Can we compare the observed periodicity with the theoretical ones? 

In order to answer these questions, Section[2]briefly reviews the existing knowl- 
edge of the chord's length in the presence of a given distribution of spheres. 
Section [3] derives two new equations for the chord's distribution in a PVT or 
NPVT environment. Section|4]is devoted to a test of the new formulas against 
a real astronomical slice and against a simulated slice. 



2. THE BASIC EQUATIONS 

This Section reviews a first example of deducing the average chord of 
spheres having the same diameter and the general formula for the chord's 
length when a given distribution for the spheres' diameters is given. 
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Fig. 2. The section of an intersected sphere of unit radius. The chord is drawn 
with the thicker line and marked with I. 
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2.1. The simplest example 
We calculate the average length of all chords of spheres having the same 
radius R. Figure [2] reports a section of a sphere having radius R and chord 
length I. The Pythagorean theorem gives 

I = ^2(1 + y)R , (1) 

and the average chord length is 

< l >=\J + i V2(l + y)Rdy = |iZ . (2) 

2.2. The probabilistic approach 
The starting point is a probability density function (PDF) for the diameter 
of the voids, F(x), where x indicates the diameter. The probability, G(x)dx, 
that a sphere having diameter between x and x + dx intersects a random line 
is proportional to their cross section 

G(x)dx l* 2 n*)dx x 2 F(x)dx 

G(x)dx ^ lx2p{x)dx ~ <x2> ■ ( 3 ) 

Given a line which intersects a sphere of diameter x, the probability that the 
distance from the center lies in the range r,r + dr is 

2nrdr 

Kr) = ^-r . ( 4 ) 

4 

and the chord length is 

I = Vx 2 - 4r 2 , (5) 

see Figure [3| 

The probability that spheres in the range (x, x + dx) are intersected to 
produce chords with lengths in the range (1,1 + dl) is 

G(x)dx 2l ^ = ^ l ^F(x)dx . (6) 
x z < x z > 

The probability of having a chord with length between (I, I + dl) is 

21 f°° 

9(1) = — / F(x)dx . (7) 
< x > Ji 

This integral will be called fundamental and the previous demonstration has 
been adapted from Ruan et al. ( 1988 ). A first test of the previous integral can 
be done inserting as a distribution for the diameters a Dirac delta function 

F(x) = 5(x-2R) . (8) 

As a consequence, the following PDF for chords is obtained: 

9(0 = ^ , (9) 

which has an average value 

<i>=^R ■ (io) 

We have therefore obtained, in the framework of the probabilistic approach, 
the same result deduced with elementary methods, see Equation 
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3. VORONOI DIAGRAMS 
This section reviews the distribution of spheres which approximates the 
volume distribution for PVT and NPVT, explains how to generate NPVT 
seeds, and derives two new formulas for the distributions of the chords. 

3.1. PVT volume distribution 



We analyze the gamma variate H(x;c) (Kiang (1966)) 

H(x;c) = -^-(cx) c - 1 exp(-cx), (11) 

r( c ) 

where < x < oo, c > 0, and T is the gamma function. The Kiang PDF has 
a mean of 

M =l, (12) 



and variance 



„> = i. (13) 



A new PDF due to Ferenc & Neda (2007) models the normalized area/volume 
in 2D/3D PVT 



3d-l 



FN(x; d) = C x x — exp (-(3d + l)x/2), (14) 
where C is a constant, 

0= , d5) 

2 2 3 / 2d (3d + ir 3/2d r(3/2d+l/2) 

and d(d = 1, 2, 3) is the dimension of the space under consideration. We will 
call this function the Ferenc-Neda PDF; it has a mean of 

M = l, (16) 

and variance 

°* = 53TT < 17 » 

The Ferenc-Neda PDF can be obtained from the Kiang function ( |Kiang 
(1966)) by the transformation 

3d + l 



2 

and as an example d — 3 means c = 5. 



(18) 



3.2. NPVT volume distribution 
The most used seeds which produce the tessellation are the so called Pois- 
sonian seeds. In this, the most explored case, the volumes are modeled in 3D 
by a Kiang function, Equation (11 1, with c = 5. An increase of the value of c 
of the Kiang function produces more ordered structures and a decrease, less 
ordered structures. A careful analysis of the distribution in effective radius 



of the Sloan Digital Sky Survey (SDSS) DR7 indicates c « 2, see Zaninetti 



(2012). Therefore the normalized distribution in volumes which models the 
voids between galaxies is 

H(x;2) = Axe~ 2x . (19) 
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3.3. NPVT seeds 

The 3D set of seeds which generate a distribution in volumes with c ~ 2 



for the Kiang function ( 19 ) is produced with the following algorithm. A given 
number, N H , of forbidden spheres having radius Rh are generated in a 3D box 
having side L. Random seeds are produced on the three spatial coordinates; 
those that fall inside the forbidden spheres are rejected. The volume forbidden 
to the seeds occupies the following fraction, /, of the total available volume 

/= 3 L3 " • (20) 

The value of c w 2 for the Kiang function is found by increasing progressively 
Nh and Rh- 



3.4. PVT chords 

The distribution in volumes of the 3D PVT can be modeled by the modern 



Ferenc-Neda PDF (14) by inserting d = 3. The corresponding distribution in 



diameters can be found using the following substitution 

x = \<\? > (21) 

where y represents the diameter of the volumes modeled by the spheres. 
Therefore the PVT distribution in diameters is 

F ^ = ^- S '"^-"""' ■ < 22 » 

Figure [4] displays the PDF of the diameters already obtained. 

We are now ready to insert in the fundamental Equation ^ for the chord 
length a PDF for the diameters. The resulting integral is 

,„ 243 Z5 2 / 3 ^ 2 / 3 81 Z 4 5 2 / 3 W/3 
9(1) 



!540 r (2/3) (e-' 3 ) 5/6 616 T (2/3) (e^ 3 ) 5/6 
135 l 7 5 2 / 3 ^6w 8 / 3 75 W/ 3 ^ 11 / 3 
+ 2464 r(2/3) (e-< 3 ) 5/6 + 4928 r (2/3) (e-' 3 ) 5/6 

125 Z 13 5 2 / 3 ^ 14 / 3 



39424 



r(2/3) (e-' 3 ) 5/6 



(23) 



This first result should be corrected due to the fact that the PDF of the 

diameters, see Figure |4j is « up to a diameter of 0.563. We therefore 
introduce the following translation 

U = L + a , (24) 



where L is the random chord and a the amount of the translation. The new 
translated PDF, gi(u;a), takes values u in the interval [—a, (6 — a)]. Due to 
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Fig. 4. PDF F(y) for the PVT diameters. 



the fact that only positive chords are defined in the interval [0, 
constant of normalization should be introduced 

<72(m; a) = Cg\(u; a), where 
C 



Jo a gi(u;a)du 



The last transformation of scale is 



R = 



and the definitive PDF for chords is 



U 



93(r;a,b) 



b 



■a)], a new 

(25) 
(26) 

(27) 
(28) 



The resulting distribution function will be 



DF lt3 (r : a,b) = g 3 (r;a,b)dr 
Jo 



(29) 



We are now ready for a comparison with the distribution function ,Fl 1 3 , for 



chord length , Li i3 in 1^,(1,3), see formula (5.7.6) and Table 5.7.4 in Okabe 



et al.| ( |2000[ ). Table [T] shows the average diameter, variance, skewness, and 
kurtosis of the already derived g 3 (r; a, b). The parameter b should match the 



average value of the PDF in Table 5.7.4 of |Okabe et aL] ( |2000[ >. The behavior 
of g${r; a, b) is shown in Figure [5j 

The behavior of DF\^(r : a, b) is shown in Figure [6j Consider a 3D 
PVT and suppose it intersects a randomly oriented line 7: the theoretical 
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TABLE 1 



THE PARAMETERS OF g 3 (r;a,b), EQ. (|28j) , RELATIVE TO THE PVT 
CASE WHEN b = 1.624 AND a = 0.563. 



Parameter 


value 


Mean 


0.662 


Variance 


0.144 


Skewness 


0.274 


Kurtosis 


-0.581 




Fig. 5. PDF gz(r; a, b) (full line) for chord length as a function of r when b = 2.452, 
a — 0.563, which means < r >= 1, and the mathematical PDF (dashed line) as 
extracted from Table 5.7.4 in Okabe et al. (2000), PVT case. 
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Fig. 6. Distribution function DF\^(r : a, b) (full line) for chord length as a function 
of r when b — 2.452, a = 0.563, which means < r >= 1, and the mathematical DF 
(dashed line) as extracted from Table 5.7.4 in Okabe et al. (2000), PVT case. 



distribution function DF\^(r : a, b) as given by Eq. (29) and the results of a 



numerical simulation can be compared. We start from 900000 3D PVT cells 
and we process 9947 chords which were obtained by adding together the results 
of 40 triples of mutually perpendicular lines. The numerical distribution of 
Voronoi lines is shown in Figure [7] with the display of DF\^{r : a,b). 



3.5. NPVT chords 



In the case here analyzed, see the PDF in volumes as given by Equa- 



tion(|19|), the NPVT distribution in diameters which corresponds to c = 2 in 
F(y) = l/3n 2 y 5 e- 1/37:y3 , (30) 



volumes is 



where y represents the diameter of the volumes modeled by spheres. We now 
insert in the fundamental Equation Q for the chord length a PDF for the 
diameters as given by the previous equation. The resulting integral, which 
models the NPVT chords, is 



gNPvr(l) = 



/7r 8 / 3 ^3(3 + 7r/ 3 ) 



5r(2/3)7r 2 ^e^ 
We now make the following translation 



(31) 



U = L + a 



(32) 




X 



Fig. 7. Comparison between data (empty circles) and theoretical curve , DFi t s(r '■ 
a,b), (continuous line) of the chord length distribution when b = 2.452, a — 0.563, 
which means < r > —1; PVT case. The maximum distance between the two curves 
is d max = 0.01. 



The new translated PDF is gNPVT,i(u; a) and the new constant of normaliza- 
tion is 

9NPVT,i{ u '-> a ) = CNPVT9NPVT,i(u;a), where (33) 

Cnpvt = -r^a : — (34) 

J gNPVT,i(u;a)du 



The transformation of scale is 



R=j , (35) 
and the definitive PDF for NPVT chords is 

, ,s 9 n pvt, 2ij, ;«) ,o R ^ 
gNPVT,3{r;a,b) = . (36) 

The behavior of the NPVT g N pvT,z{r; a, b) is shown in Figure IsJ Table [2] 
shows the average diameter, variance, skewness, and kurtosis of the already 
derived NPVT 9NPVT^{ r 'i b). The resulting distribution function with scale 
will be r 

DF NPV T(r : a,b) = / g NPVT . 3 (r;a,b)dr . (37) 
Jo 

Also in this case we produce 900000 3D NPVT cells and we process 9947 
chords which were obtained by adding together the results of 40 triples of 
mutually perpendicular lines. The numerical distribution of Voronoi lines is 
shown in Figure [9] with the display of DFnpyt(t : a, b). 
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Fig. 8. PVT PDF gs(r;a,b) for chord length as a function of r, when b — 2.452, 
a = 0.563 (full line), and NPVT PDF g NP vT,s(r; a , b ) when b = 1-741, a = 0.36 
(dashed line). In both cases < r >=1. 



TABLE 2 



THE PARAMETERS OF g nJwt , 3 (r; a, b), EQ. @, RELATIVE TO THE 
NPVT CASE WHEN b = 1.153, a = 0.36. 



Parameter 


value 


Mean 


0.662 


Variance 


1.153 



Skewness 0.324 



Kurtosis -0.442 
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Fig. 9. Comparison between data (empty circles) and theoretical curve (continuous 
line) of the chord length distribution when b = 1.741, a — 0.36, which means < r > 
=1; NPVT case. The maximum distance between the two curves is d m ax = 0.03. 
The fraction of volume forbidden to the NPVT seeds is / = 16%. 



4. ASTROPHYSICAL APPLICATIONS 

At the moment of writing, it is not easy to check the already found PDFs 
for the chord distribution in the local Universe. This is because research has 
been focused on the intersection of the maximal sphere of voids between galax- 
ies with the mid-plane of a slab of galaxies, as an example see Figure 1 in Pan 
et al. (2012). Conversely, the organization of the observed patterns in slices 



of galaxies in irregular pentagons or hexagons having the property of a tessel- 
lation has not yet been developed. We briefly recall that in order to compute 
the length in a given direction of a slice of galaxies the boundary between a 
region and another region should be clearly computed in a digital way. In 
order to find the value of scaling, b, which models the astrophysical chords, 
some first approximate methods will be suggested. A first method starts from 
the average chord for monodispersed bubble size distribution (BSD) which are 
bubbles of constant radius i?, see ^ , and approximates it with 

4 

' • (38) 



< I >= 



<R> 



where R has been replaced by < R >. We continue by inserting < R >= 
18.23/i -1 Mpc, which is the effective radius in SDSS DR7, see Table 6 in in 
Zaninetti (2012), and in this case b = 42.32/i -1 Mpc. A comparison should be 



made with the scaling of the probability of obtaining a cross-section of radius 
r, which is 31 , 33 Mpc, see Zaninetti (2012). The PDF of the chord conversely 



has average value given by the previous equation (38) and therefore b is 4/3 
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200 

Mpc/h 



400 



Fig. 10. Cone diagram of ESP galaxies when the X and Y axes are expressed 
in Mpc and H — 100 km s~ Mpc -1 (the Hubble constant), which means h=l. 
Two lines in different directions are drawn with the lengths of the chords randomly 
generated with NPVT PDF g N pvT,3(r; a, b) when b = 42.32/iT 1 Mpc and a = 0.36. 
Each chord has a different color, red and green, and different line styles, full and 
dotted. 



bigger. 

We now report a pattern of NPVT series of chords superposed on a astro- 
nomical slice and a simulation of an astronomical slice in the PVT approxi- 
mation. 

4.1. The ESP 

As an example, we fix attention on the Eso Slice Project (ESP) which 



covers a strip of 22(RA) x 1(DE) square degrees, see Vettolani et al. (19981. 
On the ESP we report two lines over which a random distribution of chords 
follows the NPVT equation ^ with b = 42.32ft,- 1 Mpc and a = 0.36, see 
Figure [10] 

4.2. A simulated slice 

In order to simulate the 2dF Galaxy Redshift Survey (2dFGRS), we sim- 
ulate a 2D cut on a 3D PVT network organized in two strips about 75° long, 
see Figure [TT] 
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Fig. 11. Portion of the NPVT V p (2, 3); cut on the X-Y plane when two strips of 75° 
are considered. The parameters of the simulation, see Zaninetti (20101, are pixels 
= 1500, side = 131 908 Km/sec and amplify = 1.2. Along a line which crosses the 
center, the chords between a face and the next one are drawn. Each chord has a 
different color, red and green, and different line styles, full and dotted. The fraction 
of volume forbidden to the NPVT seeds is / = 16% and the number of the seeds, 
77 730, is chosen in order to have b — 42.32fo -1 Mpc. 
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Fig. 12. Histogram of the frequencies in chord lengths along 31 lines (1617 sequen- 
tial chords) with a superposition of the theoretical NPVT PDF as represented by 
Eq. pit. The number of bins is 30, the reduced \ 2 is 9.21, and b = 42.32/i _1 Mpc. 



Figure [12] shows a superposition of the numerical frequencies in chord 
lengths of a 3D PVT simulation with the curve of the theoretical PDF of 
PVT chords, g(l, b), as given by Eq. (28). The distribution of the galaxies as 
given by NPVT is reported in Figure [l3j conversely Figure 15 in Zaninetti| 
(I2010I) was produced with PVT. 



5. CONCLUSIONS 

A mathematical method was developed for the chord length distribution 
obtained from the sphere diameter distribution function which approximates 
the volume of the PVT. The new chord length distribution as represented 



by the PVT formula (28) can give mathematical support for the periodicity 
along a line, or a cone, characterized by a small solid angle. At the same time, 
a previous analysis has shown that the best fit to the effective radius of the 
cosmic voids as deduced from the catalog SDSS R7 is represented by a Kiang 
function with c « 2. Also in this case we found an expression for the chord 



length given by the 3D NPVT formula (31). In order to produce this kind 
of NPVT volumes, a new type of seed has been introduced, see Section 3.3 
These new seeds are also used to simulate the intersection between a plane and 
the NPVT network. A careful choice of the number of seeds allows matching 
the simulated value of the scaling b with the desired value of b = 42.32ft -1 
Mpc, see Figure [T2| 



STATISTICS OF THE VOIDS 



17 



i — I — i — | — i — • — i — i — | — i — ■ — • — I — | — t — I — i — r 




I I I — . — ■ I ■ — ■ ■ i — I — I I i — i — I — • — i — I I I 

-0.1 o ai 

redshift 

Fig. 13. Polar plot of the pixels belonging to a slice 75° long and 3° wide generated 
by NPVT seeds. This plot contains 40 000 galaxies, the maximum frequency of 
theoretical galaxies is at z = 0.056. In this plot, Mq = 5.33 and h=l. 
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